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In this paper we investigate an approach for the numerical solution of differential equa- 
tions which is based on the perfect discretization of actions. Such perfect discretizations 
show up at the fixed points of renormalization group transformations. This technique 
of integrating out the high momentum degrees of freedom with a path integral has 
been mainly used in lattice field theory, therefore our study of its application to PDE's 
explores new possibilities. We calculate the perfect discretized Laplace operator for sev- 
eral non-trivial boundary conditions analytically and numerically. Then we construct a 
parametrization of the perfect Laplace operator and we show that with this parametriza- 
tion discretization errors - or computation time - can be reduced dramatically compared 
to the standard discretization. 

Keywords: Partial Differential Equations, Boundary Value Problems, Renormalization 
Group Methods 

1. Introduction 

It is a standard procedure to introduce a discretized space-time in the numerical 
study of physical problems. This is the case in quantum field theories where the 
numerical simulations are performed on a space-time lattice, or in classical field 
theories (hydrodynamics, electrodynamics etc.), where the corresponding partial 
differential equations can be discretized. The underlying mesh has a finite lattice 
spacing a introducing a discretization error whose size depends on a/£, where £ is 
a typical length scale of the problem. This systematical error is influenced by the 
discretized form of the differential operators. For example, the standard nearest- 
neighbor discretization of the Laplace operator has an 0(a 2 ) error. By adding 
additional couplings this error can be reduced to C(a 4 ), which makes it easier to 
perform the continuum limit a — > 0. In principle, the discretization error can be 
eliminated order by order this way, but the resulting Laplace operator would not 
have much practical value in general: it would become a broad, non-local operator. 
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The notion of locality will play an important role in the following. We shall 
call a discretized differential operator local if the size of the region R where the 
couplings are significantly different from zero goes to zero relative to the typical 
length scale £ of the problem in the continuum limit: — ► (a — > 0). The 

standard nearest-neighbor discretization of the Laplace operator is obviously local. 
Similarly, a lattice difference operator whose couplings decay exponentially oc e ,r 
with the distance r between the connected points, where 7a = 0(1), is local also. 
In this case oc 1/(7^) oc a/£ — * in the continuum limit. 

Is it possible to construct "perfect" discretized differential operators, i.e. opera- 
tors which are local and lead to differential equations on the lattice whose solution 
is free of discretization errors even if a/£ is not small and thus the resolution is 
bad? The answer is yes. The problem is related to the existence of perfect lattice 
actions of classical field theories 111. The Euler-Lagrange equations corresponding to 
such actions are perfect difference equations in the sense discussed above. Their 
existence and the technical way of constructing them follows directly from Wilson's 
renormalization group (RG) theory □. Perfect classical actions were constructed 
and tested not only for free field theories like_the scalar field H, but also foi_such 
non-trivial cases like thenon-linear er-model era, Yang-Mills gauge theory ErEI or 
the Schwinger model EHj'E2l. These results might open new paths in the numerical 
study of partial differential equations. The first attempt in this direction has been 
made by Katz and Wiese for fluid dynamics til. Although the latter problem seems 
to be simpler than those treated before, there are new difficulties here — among 
them the influence of non-trivial boundary conditions. Consider the following sim- 
ple example: the potential energy of a massive membrane fixed to a frame at its 
boundary. This is a good problem for testing, as it can be solved analytically for a 
square frame of size L x L and therefore the exact solution is known. The standard 
procedure to solve this problem numerically is to define a lattice and calculate the 
amplitude of the membrane at the lattice points using the standard discretization 
of the Laplace operator. However, if a/L is not very small, the results will be far 
from the exact value due to the large discretization error. To get acceptable results, 
one needs to go to very fine lattices, that is to very high resolution. On the other 
hand, as we mentioned, the perfect lattice Laplace operator gives the exact value for 
the quantity we are looking for at any lattice size, at any resolution, even when the 
lattice consists of only one point. For the boundary conditions mentioned above, 
the perfect Laplace operator can even be found analytically. 

So far, this is theoretical. For practical calculations, one should use a trun- 
cated operator which is easy to handle. For example, one might consider only 
nearest-neighbor and next-to-nearest-neighbor couplings between lattice points. A 
boundary will then influence these couplings. Such a truncation only works when 
the neglected couplings are small, that is when the operator is local. If the influence 
of the boundary is also local, a parametrization could be constructed which can be 
used for arbitrary boundary shapes. In this paper, we construct and test such a 
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parametrization. Let us summarize the problems addressed: 

• How do non-trivial boundary conditions affect the perfect Laplace operator? 

• Can we truncate and parametrize a perfect difference operator in a way that 
it is useful for practical purposes? 

• Is the improvement we achieve for specific results when using a parametrized 
perfect operator worth the effort? 

As this is a field practically untouched, we turn to the basic problems first. We will 
see that there the concept of perfect actions leads to excellent results. 

2. Renormalization Group Transformations and Perfect Actions 

Let us first give a brief review of the construction of perfect actions. Consider a 
theory described by an action A((j) n ), where <j) n is the (i-dimensional field variable 
at the lattice point n — (m, . . . , rid). On this lattice, form blocks of 2 d points each. 
For every block, define a blocked field variable 

nen B 

If the original lattice spacing has been a, the blocked field Xn B lives on a lattice 
with spacing 2a. A Renormalization Group transformation (RGT) step leads to 
a new action A'(x) for the blocked variable by integration over the original field 
variables (f> n : 

e- A ' M =Uf d^Yl d (Xn B -by d J2 Me~ AW - (2.2) 

It remains to rescale the unit of length in order to keep physical length scales 
unchanged. After such an RGT step, our theory still describes the same long- 
distance behaviour, i.e. the same physics. But we have reduced the number of 
degrees of freedom, that is the number of space-time variables, by a factor of 2 d . On 
the other hand, the form of the action has changed. Iterating this transformation, 
one gets coarser and coarser lattices without adding new discretization errors. 

An interesting property of RG transformations is the occurrence of fixed points. 
Most generally, the action of a theory consists of all kinds of interactions and can 
be written as a sum of interaction terms 0i(<j>), i = 1, 2, . . . 

A{<f>)=Y / K i 6 i {<t>), (2.3) 

i 

where Ki are the respective coupling constants. As we said, a RG transformation 
changes the form of the action, so repeated transformations generate a flow in 
coupling constant space 

{ K W}^{kW}^{kW}^... . (2.4) 



4 S. Hauswirth 



A fixed point occurs if the set of coupling constants remains unchanged under a 
RG transformation {K^} = {JQ n = {A",*}. The fixed point action, which we 
designate by an asterisk, is then defined as 

^(0)=^A*^(0), (2.5) 

i 

and depends on the explicit form of the RG transformation. The fixed point ac- 
tion has the beautiful property of being classically perfect: It reproduces all the 
important classical properties of the continuum action EJ. 

Now let us consider an example: The perfect Laplace operator A* on a d- 
dimensional space-time without boundaries can be calculated from the fixed point 
action of a free real scalar field with the continuum action 

-4(0) = \j ddx d^{x)d^{x). (2.6) 

The equation of motion for this action is the Laplace equation. A general discretiza- 
tion of the action (|j]) contains terms which couple the field at one lattice site to 
the field at another one: 

A{(f>) = - J2<l>np{r)<f>n+r, (2.7) 
n,r 

with the coupling constants p(r), r — (n, ... , frf). For the standard Laplacian in 
two dimensions, p(0) = 4 and p(r) = —1 (V|r| = 1). We generalize the block 



transformation in (2^) and write 

^*)=n fd4 n e- A ^- T ^\ (2.8) 



c • e 



where c is a normalization constant and T(x, <f) is the blocking kernel 

n B nGriB 

The parameter n gives us the possibility to optimize the fixed point action for 
locality. For k — > oo, this RG transformation goes over tQihe one in Eq. (2.2). The 



to the 

□d. 



fixed point of Eq. ( J2.8| ) can be calculated analytically EjO. In momentum space, 
the inverse of p* (q) is 

1 v ! TT sin2(f + 7rU i 1 (210) 



The fixed point Laplacian (A*0)„ = — ^ r p*(r)(f> n+r is then calculated numerically 
by Fourier transforming p*(q). A thorough examination of different blocking kernels 







has shown that the kernel (2.9) with k — 2 gives very good results in terms 
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of locality: The couplings p(r) decay exponentially oc e~ 7r with a large decay 
coefficient 7 w 3.5. For d = 2, we get the translationally invariant perfect Laplacian 
on a plane, which we will use in the next section. 

3. Perfect Laplacian for a Square Boundary 

The perfect Laplace operator calculated above only holds in the absence of bound- 
aries. So let us now introduce non-trivial boundary conditions: Restrict the field <f> 
to a square area in two dimensions with (/) = on the boundaries. In the presence of 
boundaries, the perfect Laplacian p*(n,n') will no longer be translation invariant. 
In this section we present two ways of calculating the perfect Laplacian for these 
boundary conditions. First we explicitly perform the RG transformation, leading 
to an analytical expression for the perfect Laplacian, and second we show how to 
find p*(n,n') as a function of the translationally invariant perfect Laplacian p*(r), 



which is already known from ( 2.10 ). The latter procedure will lead to an elegant 



way of avoiding possible normalization problems when truncating the couplings. 
3.1. Explicit RG Transformation 

Consider a square box with side length L in d = 2. In the continuum, the action of 
a free scalar field is 

A<t>) = l [ d 2 x d^{x)d^{x). (2.11) 
1 Jo 

For the calculation of the fixed point action, we will work in Fourier space. Our 
ansatz is to use a Fourier transformation 

0(*) = 7aZ)*9( a: )te ( 2 - 12 ) 



with orthonormal basis functions & q (x) = 2 sin(gixi) sin^a^) which ensure that 
the field <f> fulfills the boundary conditions. The momentum variable q = (91,92) 
takes the discrete set of values = kiir/L with fej = 1,2,..., 00 (i = 1,2). In 



momentum space, the action (2.11) reads 

As shown in Fig. [l], we define the lattice points n, at half-integer values of the 
continuum space variable x% = (m + l/2)a with < Xi < L = Na. Performing the 
RG transformation ( |2.8| ), ( [2.9| ) on this lattice, we arrive at a coarser lattice with 
lattice unit 2a, defined in the middle of the dashed blocks in Fig. [|. Consider the 
correlation function 



# n e-^V n ,0 n „, (2.14) 
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Figu re 1: The TV x TV-lattice in an area with square boundaries. The blocking kernel 
(2.9) connects the four fine lattice points inside a dashed box to one coarse lattice 
point in the center of the box. 



where Z is the partition function and assume A = A*. From Eq. (2.8) follows that 
the correlation functions on the coarse and fine lattices are related by 



(Xn B Xn' B ) = (j) £ £ Mn-) + ^ 



(2.15) 



h£hb n'Gn' 



Iterating this transformation an infinite number of times leads to a relation between 
the correlation functions on the original lattice with lattice unit a and the one in 
the continuum. (As we need no longer two lattice field variables, we replace the 
coarse lattice variable \n B by our standard lattice notation <f> n again.) 



l r i 

2„ / jl j I 



1 



{<t>n4>n') = d x Id x'{4>((n + x)a)4>((ri + x')a)) + — <W- 
Jo Jo 3k 

On the right hand side, we insert the continuum free field propagator 



(2.16) 



(2.17) 



which follows from (2.13), with x^yi G (0, L). Setting the lattice unit a — 1 and 
performing the definite integration over x and x' gives 



4 2 1 1 2 (sin 22 -) 2 1 1 

(0n<M = 2J II sin (li( n i+2^ Sm faK + ^)) II (21)2 2+^' 5 " 

After a few algebraic steps, we can bring this into the following form: 



1 1 



n 



(2.18) 



2 sin 2 Si 1 

3; , 2 „ +— <W, (2.19) 



(g + 27r0 2 ^ (§1+^)2 3k 
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with ^ Q (n) = ^(nx) ■ £,Q 2 (n 2 ), and 

e C« ^ - / V2sm(Qi(rii + |)) ,Qi^n, ,„ 

sinCQ^ + i)) , Qi = 7t. ^- 2UJ 

The new momentum variable Q is restricted to the Brillouin zone and takes the 
values Qi — kiir/N, ki = l,..,,N. The result (|2.19| ) provides us the fixed point 
propagator in momentum space 

1 00 1 2 .2 0; 1 

1 x - 1 tt sm ^ 1 



x ^ 1 -1— 1- aiu. — 1 

Perform the Fourier transformation 

*') = ^2 E ' ( 2 - 22 ) 



to get the couplings of the fixed point action A*{4>) — 1/2 , p*(n, n' 
in configuration space (rij = 0, 1, . . . , AT — 1). In contrast to the situation on an 
unbounded lattice, p*(n, n') is no longer translationally invariant; it does not depend 



only on the distance r = n' — n as the couplings in Eq. (2.7). Consequently, we 
have a different p*(n,r) for every lattice point n. On the other hand we expect 
— and it is really the case, as we will show — that away from the boundaries the 
difference between p*{n, r) and the perfect Laplacian of the translationally invariant 
case p*(r) goes to zero exponentially with the distance from the boundary, which 
will be a crucial property to construct a parametrization. 

3.2. Construction from Symmetry Properties 

There is another, more elegant way to find the perfect Laplacian for a square region 
with zero boundary conditions in d = 2. The translationally invariant perfect 
Laplacian p*(r) ( [2.1C| ) is already known, and we show that this knowledge can 
be used to solve the problem with these non-trivial boundary conditions. As a 
boundary, consider a single wall at xi = first where the field <fi has to vanish. We 
extend the field beyond the boundary and introduce the condition 

4>(-Xi,x 2 ) = -<j>{xi,X2), (2.23) 

which implies <f> = on the boundary. For the lattice field, we have 0_ ni _i !ri2 = 
— (f> ni ,n 2 y as the lattice points n are at half-integer values of x = n+ 1/2. Using this 
symmetry relation and the symmetry property 7*1, 7*2) = p*{i~i, r 2 ), we rewrite 
the perfect Laplace equation on the unbounded lattice into an equation on the right 
halfplane: 

]Tp*(r)0 Il+r - ^[^p*(r-n)0 r +^p*(r-n)0 r ] (2.24) 

r r 2 n>0 ri<0 

= ^2 ^2 [p*(r -n) - p*(r a +ni + l,r 2 - n 2 )]<j) r . 

r2 ri>0 
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1 1 


Pi(n,r2) 


n = -2 


n = —l 


n = 


n = 1 


n = 2 


;. . . .; 


r 2 = 2 


0.00162 


-0.00068 


-0.00200 


-0.00068 


0.00162 


i i 


T2 = 1 


-0.00068 


-0.19024 


-0.61773 


-0.19024 


-0.00068 


pi 


r 2 = 


-0.00200 


-0.61773 


3.23881 


-0.61773 


-0.00200 




r 2 = -1 


-0.00068 


-0.19024 


-0.61773 


-0.19024 


-0.00068 




r 2 = -2 


0.00162 


-0.00068 


-0.00200 


-0.00068 


0.00162 


















P2 (1-1,1-2) 


n = -2 


n = -1 


n = 


n = 1 


ri = 2 


!: : : : :! 


r 2 = 2 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


■ i 


r 2 = 1 


-0.00230 


-0.18956 


-0.61573 


-0.18956 


-0.00230 


P2 


r 2 = 


-0.00200 


-0.61773 


3.23881 


-0.61773 


-0.00200 




1"2 = —1 


-0.00068 


-0.19024 


-0.61773 


-0.19024 


-0.00068 




r 2 = -2 


0.00162 


-0.00068 


-0.00200 


-0.00068 


0.00162 



P3(n,r2) 


ri = -2 


ri = -1 


n =0 


n = l 


ri = 2 


r 2 = 2 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


r2 = 1 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


r 2 = 


-0.00132 


-0.42749 


3.85654 


-0.42749 


-0.00132 


1"2 = —1 


-0.00230 


-0.18956 


-0.61573 


-0.18956 


-0.00230 


r 2 = -2 


0.00162 


-0.00068 


-0.00200 


-0.00068 


0.00162 



P4(r i,r2) 


n = -2 


n = -1 


n = 


n = 1 


n = 2 


r 2 = 2 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


r2 = 1 


-0.00230 


-0.18956 


-0.61573 


-0.18726 


0.00000 


r 2 = 


-0.00200 


-0.61773 


3.23881 


-0.61573 


0.00000 


1"2 = — 1 


-0.00068 


-0.19024 


-0.61773 


-0.18956 


0.00000 


r 2 = -2 


0.00162 


-0.00068 


-0.00200 


-0.00230 


0.00000 



P5(n,r2) 


n = -2 


n = -1 


n = 


n = 1 


n = 2 


r 2 = 2 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


T2 = 1 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


r 2 = 


-0.00132 


-0.42749 


3.85654 


-0.42617 


0.00000 


T2 = — 1 


-0.00230 


-0.18956 


-0.61573 


-0.18726 


0.00000 


r 2 = -2 


0.00162 


-0.00068 


-0.00200 


-0.00230 


0.00000 



pe(n,r2) 


ri = -2 


r\ = —1 


ri = 


n = l 


ri = 2 


r 2 = 2 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


r2 = 1 


0.00000 


0.00000 


0.00000 


0.00000 


0.00000 


r 2 = 


-0.00132 


-0.42617 


4.28403 


0.00000 


0.00000 


1"2 = —1 


-0.00230 


-0.18726 


-0.42617 


0.00000 


0.00000 


r 2 = -2 


0.00162 


-0.00230 


-0.00132 


0.00000 


0.00000 



Table 1: The normalized, truncated couplings of the perfect Laplacian for points 
far from boundaries (pi), for points near a wall (p2,P3) and near a convex corner 
(p4,P5,Pe)- I n the small image on the left, the lattice point n is denoted by a 
circle, and all the sites in the 2-hypercube around n which lie inside the boundary 
are shown as small dots. The table on the right shows the values of the (r 1; r 2 )- 
coupling for the lattice point n. The couplings for P2 _ P6 are made up from p\ with 
the construction explained in the text. 
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This defines a perfect Laplace operator with <fi = on the second axis, that is the 
perfect Laplacian near a wall 

p*(n, r) = p*(r) - p*{r x + 2n x + 1, r 2 ), (2.25) 

where n x is the distance from the wall. The variables n and n! = n + r are here 
restricted to the area inside the boundary which is the positive halfplane in the first 
coordinate. 

Now consider a boundary of two walls at x x = and x 2 = forming a corner. 
In the Laplace equation, the sum over the lattice points can be split into sums over 
the four quadrants of the plane. Proceeding as in the above case for a wall, we get 
the perfect Laplacian 

p*(n, r) = p*{r) - p*(r x + 2n x + 1, r 2 ) - p*(r x ,r 2 + 2n 2 + 1) (2.26) 

+P*(n + 2rn + 1, + 2n 2 + 1), 

with both n and n + r lying in the first quadrant. 

Finally, we may form a square boundary out of four walls at X\ = 0, x\ = 
N, a*2 = and x 2 — N. The field outside the boundary is formally defined by 
periodically mirroring it at the boundary with alternating sign 

4>(x 1 ,x 2 ) = (~l) kl+k ^({-l) kl x 1 + 2l 1 N,{~l) k *x 2 + 2l 2 N) 1 (2.27) 

for any fej = 0, 1 and U € Z, (i = 1, 2). The sum over the whole plane then splits 
up into sums over iV-squares 

2 1 N-l 

r i=l li£Z k i= Si=0 

where the variable r running over the lattice points is given by = (— 1) *Sj — hi + 
2kN for i — 1,2. The perfect Laplacian consists of the infinite sum of couplings 

2 1 

p* {n,r)=J2J2J2 (-l) fel+ V (ri + 2fcmi + ki + 2hN, r 2 + 2k 2 n 2 + k 2 + 2l 2 N) , 

'=1 li£Z ki=0 

(2.29) 

with both n and n + r lying inside the boundary. A check for a lattice with N = 1 
shows that the relation ( [2.29 ) with the translationally invariant couplings ( 2.1 Cl| ) 
gives the same result p*(0, 0) = 4.95513 as the explicit RG transformation (2.21), 

(EEH). 

4. Perfect Laplacian near a Concave Corner 

For boundary shapes where the fixed point action is hard or impossible to find 
analytically, RG transformations may be performed numerically. We performed 
such a numerical RGT in order to find the couplings of the fixed point Laplacian 
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/J 2 



□ □ 



Resulting coarse field couplings 



Fine 



couplings needed for input 



hi 



N 



Figure 2: The N = 32 - lattice used to calculate the couplings near a concave 
corner. The couplings between the fine lattice points (shown as dots) have to be 
fixed. The RG transformation then gives the couplings between the coarse lattice 
points (denoted by squares) inside the large dashed box, which are used as an input 
for the fine field couplings inside the small dashed box in the next iteration. 



near the concave corner of an L-shaped boundary, that is a square with a smaller 
square cut out, in d = 2 (see Fig. ||). For quadratic actions and blocking kernels, 
the RG transformation Eq. (2.8) can be written as a minimizing condition for the 
fine field 6: 



A'(x) = min [A(4>) + T (x, 4>)] + const. 



(2.30) 



To get close to a fixed point, we have to iterate this RGT step. The results of the 
previous step — which are the couplings of the resulting coarse action — are then 
used as an input for the next step, that is as a new starting guess for the fine action 
A{4>). After O(20) iterations, we find a very close approximation to the couplings 
of the fixed point action. 

We worked on a lattice with N — 32, as shown in Fig. |[ and looked for the 
couplings of the coarse field inside the large dashed box. Outside the small dashed 
box in Fig. |^, we used in every RGT step the parametrized couplings from Table 
|l| for the fine action. For the couplings inside the small dashed box, we used the 
standard Laplacian in the first iteration, and afterwards the result of the previous 
iteration. For every step, the coarse field couplings can be read out one by one by 
choosing appropriate coarse field configurations as an input. Our results for the 
fixed point couplings are listed in Table @. 



5. Parametrization for Boundaries of Arbitrary Shape 
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P7 



P7(ri,r 2 ) 


n = -2 


n = -1 


n = o 


n = l 


n = 2 


r 2 = 2 


0.00160 


-0.00072 


-0.00192 


-0.00206 


0.00000 


r 2 = 1 


-0.00071 


-0.19033 


-0.61787 


-0.19145 


-0.00206 


r 2 = 


-0.00190 


-0.61784 


3.26307 


-0.61787 


-0.00192 


r-2 = —1 


-0.00071 


-0.19033 


-0.61784 


-0.19033 


-0.00072 


r 2 = -2 


0.00160 


-0.00071 


-0.00190 


-0.00071 


0.00160 



PS 





ri = -2 


n = -1 


n =0 


n = l 


n = 2 


r 2 = 2 


0.00160 


-0.00071 


-0.00043 


0.00000 


0.00000 


r 2 = 1 


-0.00071 


-0.19032 


-0.61625 


-0.18956 


-0.00233 


r 2 = 


-0.00190 


-0.61787 


3.26192 


-0.61787 


-0.00196 


r 2 = —1 


-0.00071 


-0.19033 


-0.61784 


-0.19033 


-0.00072 


r 2 = -2 


0.00160 


-0.00071 


-0.00190 


-0.00071 


0.00160 



■ ■©■ 



P9 



P9(ri,V2) 


ri = -2 


ri = —1 


ri = 


n = l 


ri = 2 


r 2 = 2 


0.00160 


0.00062 


0.00000 


0.00000 


0.00000 


r 2 = 1 


-0.00071 


-0.18948 


-0.61637 


-0.18963 


-0.00233 


r 2 = 


-0.00192 


-0.61787 


3.27774 


-0.61788 


-0.00196 


r 2 = -1 


-0.00072 


-0.19033 


-0.61784 


-0.19033 


-0.00072 


r 2 = -2 


0.00160 


-0.00071 


-0.00190 


-0.00071 


0.00160 



PlO 



pio(ri,r 2 ) 


ri = -2 


n = — l 


ri = 


ri = 1 


ri = 2 


r 2 = 2 


0.00021 


0.00000 


0.00000 


0.00000 


0.00000 


r 2 = 1 


-0.00206 


-0.18964 


-0.61594 


-0.18965 


-0.00233 


r 2 = 


-0.00196 


-0.61788 


3.24479 


-0.61788 


-0.00196 


r 2 = — 1 


-0.00072 


-0.19033 


-0.61784 


-0.19033 


-0.00072 


r 2 = -2 


0.00160 


-0.00071 


-0.00190 


-0.00071 


0.00160 



Pn(Ti,r 2 ) 


n = -2 


n = — l 


n = o 


n = l 


n = 2 


r 2 = 2 


0.00155 


-0.00206 


-0.00051 


0.00000 


0.00000 


r 2 = 1 


-0.00068 


-0.18948 


-0.46147 


0.00000 


0.00000 


r 2 = 


-0.00186 


-0.61625 


3.36152 


-0.46147 


-0.00051 


r 2 = — 1 


-0.00076 


-0.19145 


-0.61625 


-0.18948 


-0.00206 


r 2 = -2 


0.00161 


-0.00076 


-0.00186 


-0.00068 


0.00155 



P12 



Pi 2 (»*i,r 2 ) 


n = -2 


ri = -1 


n = o 


n = l 


n =2 


r 2 = 2 


0.00021 


0.00055 


0.00000 


0.00000 


0.00000 


r 2 = 1 


0.00062 


-0.02666 


0.00000 


0.00000 


0.00000 


r 2 = 


-0.00043 


-0.46147 


3.88365 


-0.42738 


-0.00121 


r 2 = — 1 


-0.00206 


-0.18956 


-0.61637 


-0.18964 


-0.00233 


r 2 = -2 


0.00156 


-0.00074 


-0.00179 


-0.00073 


0.00155 



P13 



Pia(ri,r 2 ) 


n = -2 


n = — l 


n = o 


n = l 


n = 2 


r 2 = 2 


-0.00001 


0.00000 


0.00000 


0.00000 


0.00000 


r 2 = 1 


0.00055 


0.00000 


0.00000 


0.00000 


0.00000 


r 2 = 


-0.00051 


-0.42738 


3.85640 


-0.42757 


-0.00121 


r 2 = — 1 


-0.00233 


-0.18963 


-0.61594 


-0.18965 


-0.00233 


r 2 = -2 


0.00155 


-0.00073 


-0.00179 


-0.00073 


0.00155 



Table 2: The truncated couplings of the perfect Laplacian near a concave corner. 
These couplings are the results of the iterative blocking procedure and are therefore 
not normalized. 
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The main observation which allows us to construct a reasonable parametrization for 
the perfect Laplacian that can be used to approximate arbitrarily shaped boundaries 
is the fact that not only its couplings decay exponentially with the distance, but also 



the effect of the boundary on the couplings: According to Eq. ( 2.25 ), the difference 
between the perfect Laplacian near a wall and the translationally invariant perfect 
Laplacian on a plane without boundaries 



Ap* n (r)=p*(n,r)-p*(r)=-p*(r 1 + 2n 1 + l,r 2 ), 



(2.31) 



decays with the distance from the wall n\ at twice the decay rate of the FP Laplacian 
couplings. The value of the decay constant 7, which is defined by Ap*(r) oc e~ 7 ' ni , 
where n\ denotes the orthogonal distance from the wall — or the diagonal distance 
from the corner, respectively — is listed for several couplings r near a wall and near 
a corner m Table § This very strong exponential decay means that for points a few 



-5 , \ 



o r ■ 
□ r : 
o r ■ 



(0, 1) 
(1,0) 
(1, 1) 



-10 



-20 



;25 



C 



-30 L 



4 

m 



Figure 3: The relative difference Ap* (r)/p*(r) between three couplings near a wall 
and the respective couplings of the perfect Laplacian on an infinite lattice. The 
couplings p*(n,r) are calculated on a lattice with N — 17 at the lattice sites n = 
(ni, (N + l)/2). n\ is therefore the distance from the wall in lattice units. The 
difference decreases exponentially with the distance from the wall. The slope 7^ is 
listed in Table |3[ For distances larger than 3 lattice units, the difference is already 
beyond the numerical accuracy. 

lattice spacings away from walls, the presence of the boundary influences the fixed 
point Laplacian at that point only negligibly, and therefore we may parametrize the 
perfect Laplacian in a way that it can be used for lattices with boundaries of an 
arbitrary shape. 
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r 


^(wall) 


^{corner) 


(0,1) 


7.2 


7.2 


(1,0) 


6.5 


7.2 


(1,1) 


7.8 


7.0 


(0,2) 


5.9 


6.0 



Table 3: The exponential decay rate 7 for the effect of the boundary on different 
couplings. The values are from linear least-squares fits to the first few data points 
in Figure [| 



In our parametrization, we take together points with very similar p*(n,r) (e.g. 
all the points with a distance from the boundary of more than one lattice unit) and 
use the same operator p m {r) for all these points. Thus, we classify the lattice points 
by their position relative to the boundary (lower left corner, right wall, . . . ) and 
define a set of operators p m (r), m = 1, . . . , 13, where m is the type of lattice point 
as shown in Fig. ^. Making this approximation, our set of operators is no longer 
perfect, and so we have dropped the asterisk in the notation. 



Figure 4: Point types used for the parametrization. The lattice points are classified 
by their relative position to the boundary. For every point type m we define a 
different Laplace operator p m (r) (see Table [l]). The definition of pf,...,px$ is 
shown in Table 0. 



5.1. Truncation and Normalization 

For practical purposes, the perfect Laplacian must be truncated to a finite number 
of couplings. The most primitive nearest-neighbor Laplace operator and its simple 
improved versions have the basic property that they all approach the same exact 
continuum result as the resolution is increased (universality). We have to ensure 
that our approximate truncated perfect Laplace operator also has this basic prop- 
erty. This can be achieved when imposing some conditions implied by elementary 
principles on our parametrized action. Let us consider the case of an unbounded 
lattice first. Any discretized action should describe the physical properties of the 
continuum action. For example, the spectrum of the discretized action, which is 
given by the poles of the propagator, has to go over to the spectrum of the contin- 
uum action. In the continuum, the free field propagator is q~ 2 , and the poles are 
at q = (p, i\p\)' The relativistic dispersion relation in two euclidian dimensions is 
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E(k) = —iq 2 = \p\. The propagator l/p(q) of a quadratic discretized action should 
have the same poles, therefore the couplings have to fulfill the sum rule 

S a = ^p(n,r 2 ) = 0. (2.32) 

ri,r 2 

Furthermore, we demand that the lattice action takes the form of the continuum 
action for small momenta g, that is lim^o p{q) — <7 2 - This leads to the sum rule 



(r 1 +r|)-p(ri,r 2 ) = -4. (2.33) 



When we have boundaries, the above considerations can't be used directly to 
find sum rules. But with a different approach, we can at least find some conditions 
for the couplings. Consider a wall at x\ = 0. The field 4>{x) = x\ is a solution to 
the continuum Laplace equation and is zero at the boundary. Therefore the lattice 
field <&„ = rii + 1/2 is a solution to the lattice Laplace equation, and when inserting 
this solution we find a condition for the couplings p n {r): 

S{ w) = j2p-(ri,r 2 )(n + l) = 0. (2.34) 

ri,r 2 

For a corner, <fi(x) — X\X 2 is a solution in the continuum. Hence, the lattice field 
&n = (ni + l/2)(ri2 + 1/2) solves the lattice Laplace equation, and the sum rule is 

S[ C) = E Pn(n,r 2 )(ri + \){r 2 + \) = 0. (2.35) 

ri,r 2 



While ( |2~34| ) and fl2~35[ ) are trivially fulfilled on an unbounded lattice, they make 



sense as a normalizing condition for the couplings in the presence of a boundary. 
We remark that Eq. ( [2.35 ) holds for both convex and concave corners. 



We construct our parametrization in the following manner: First we normalize 
the truncated Laplacian for inner points pi — which we take from Eq. ( [2.1C| ) — by 
setting the (0, 0)-coupling properly to get Sq — and by rescaling all couplings to 
ensure S2 = —4. Then we build the Laplace operators p 2 ,- ■■ ,pe out of the already 



corrected couplings of p\ with the help of the relations ( 2.25 ) and ( 2.26 ). The 
resulting couplings are listed in Table 0. Together with the couplings near a concave 
corner P7, . . . , P13, we have a parametrization that can be used to approximate any 
two-dimensional shape. 

6. Testing the Parametrization 

Let us first check the statement that the fixed point action is classically perfect on 
an exactly solvable problem. Consider the free scalar field with a constant source 

A<f>) = \ f d 2 x d^{x)d^{x) + f f d 2 x^{x). (2.36) 



Jo 
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We may interpret this as the action of a membrane which is fixed to <f> — at the 
boundaries and has mass density /. The total mass is then F = fL 2 . The equation 
of motion to the action ( p.36| ) is the Poisson equation Acj>(x) — f. The Green's 
function to the Laplace operator which fulfills the given boundary conditions is 

G(I , 9) = ^ E iM^>, , 2 . 37) 



with ^ q (x) — 2sin(gicci) sin(g 2 ^2) and = kiir/L, (fcj € N). Plugging ( |2.37| ) 



back into ( |2.36| ), we can calculate the value of the action and the potential energy 
E (cont) _ 2A _ _ .035144 • F 2 . For k ->• oo, it is easy to check that the perfect 



action (2.22) reproduces this continuum result exactly even for N = 1, i.e. when 
the lattice consists of only one single point! 

To compare classical quantities like the value of the action on the lattice and in 
the continuum, we have to take into account the constant c which appears in the 



RG transformation (2.8). For the fine lattice action 

^) = ^5>*(r)<M„ +r + /^5>„, (2-38) 



n , r 



with the perfect Laplace operator p*(r) and constant mass density fy, we get after 
one RGT step the coarse lattice action 



1 f 2 

■A(X) = ^ P*( r B)Xn B Xn B +r B + SxzZXtlB ~z2^ ( 2 - 39 ) 

n B ,r B 7i B n B 

with f x = 4f<p. Hence, after iterating to the continuum, the potential energy is 
calculated from the lattice field configuration <fi which solves the perfect Poisson 
equation by 

n 

where V denotes the lattice volume. 

As a check of the quality of our parametrization of the perfect Laplace operator, 
we numerically solved the lattice Poisson equation A<f> — f for square boundaries 
with the standard (4, —1, — 1, — 1)-Laplacian and with our parametrized Laplace 
operator made up from pi,...,pQ. For the resulting field configuration, we de- 
termined the potential energy E — f ' for the field configuration 0( s * d ) 



computed with the standard Laplacian and (2.4C) for the configuration 0( pp ) com- 
puted with the parametrized perfect Laplacian. Fig. ^| shows a plot of the relative 
error (E — £ , ( cont ))/i*;( co ™ t ) as a function of the inverse lattice volume 1/N 2 . The 
error for our parametrization is proportional to 1/N 2 and is for any lattice size 
smaller by a factor of about 180 than the error of the standard Laplacian. This 
factor can probably still be increased by tuning the normalization procedure of the 
truncated couplings. Fig. ^| shows what happens when the couplings are not nor- 
malized and therefore do not fulfill the sum rules So and Sa- We see that the error 
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is then no longer proportional to the inverse lattice size. While for coarse lattices 
— that is small N — the error is small, the results get worse for finer lattices, and 
somewhere a point is reached where the standard Laplacian gets better than the 
paramctrization. If this should be avoided, one has to correct the truncation error. 



0.020 
0.015 
0.010 

U 

W 

"0.005 



& - -0 Standard Laplacian 
o-^Param. Perfect Laplacian 



-ms®®&x&$$-& Q-<$- e--® ^ 







0.01 



1/N 2 



0.02 



Figure 5: Relative error for the numerical solution of our test problem. The error 
for the parametrized perfect Laplacian is for any resolution about 180 times smaller 
than for the standard Laplacian and gives a 0.06% error at N = 5. In other words, 
our paramctrization gives better results on a 5 2 lattice than the standard Laplacian 
on a 74 2 lattice. 



Finally, we make the same comparison for L-shaped boundaries as shown in 
Fig. ||. As above, we solve the lattice Poisson equation with the standard Laplacian 
and the parametrized perfect Laplacian, for which we take the couplings pi , . . . , P13 
listed in Tables | and §. The lattice volume V for this shape is V = 3/47V 2 . We 
compare lattice sizes N between 8 and 72. As we see in Fig. [?], the parametrized 
perfect Laplacian gives excellent results for any lattice size, while for the standard 
Laplacian we have to use very fine lattices to get acceptable results. For example, 
take a lattice with N = 10. The standard Laplacian then has an 0(25%) error, 
while the error for the parametrized perfect Laplacian is O(0.1%). 

7. Generalizations 

The parametrization presented in the previous sections has some limitations: The 
boundaries have fixed orientation and are always half a lattice spacing away from 
the lattice points, and we only used a constant source in the action. In this section 
we show how to overcome these limitations. 
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Figure 6: Relative errors for the correctly normalized parametrized perfect Lapla- 
cian and for an incompletely normalized parametrization. Although the correct 
normalization gives somewhat larger errors for small lattices (large lattice unit a), 
it makes sure that for any lattice size, the results are better than for the standard 
Laplacian. The error of the standard Laplacian is far beyond the scale of this plot 
(compare Figure ||) . 

Consider a RG transformation of the action with a general source term J{x): 
1 



S.M = § / ^ l d ^( x ) d ^( x ) + J ( x )0( x )]- ( 2 - 41 ) 
Blocking out of continuum 

an 

with the kernel 

T H [9>,<f>] =«^Y$„- f d d x lu(x ~ n)<f>(x)y , (2.42) 
n 

where lu(x) is an arbitrary blocking function gives the the result in Fourier space 



Am = r 



d d k - 



d d k ~ 



Q(k) ~, 



— dm p*(kM-k)- J —<f>(-k)p*(k)-^J(k)-W(J*) 



(2.43) 

W{ J 2 ) denotes a ^-independent term where the source only appears in second order 
and the inverse of p*(k) is the fixed point propagator 

1 Q(k + 2irl)u>(-k - 2-kI) 1 

(fc + 27T/) 2 + k' 



p*(k) 



iez d 



(2.44) 



The equation of motion to the fixed point action ( |2.43 ) in configuration space is the 
perfect Poisson equation 



(2.45) 
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Figure 7: Comparison of the results for the potential energy of an L-shaped mem- 
brane between standard and parametrized perfect Laplacian. The error for the 
standard Laplacian is of order 1/7V. 



with 

J FP (k)=p*(k)^-J(k), (2.46) 

Therefore the perfect Laplace operator can be used to solve the Poisson equation 
with any source J(x). 

From Eq. (2.43) we can also read the fixed point field operator (f> FP which 
remains unchanged under RG transformations. Assume that the source J(x) is 
small. Then the term of order J 2 can be neglected, and the transformation back to 
configuration space yields 



•4*1 =o / 7^mP*( k M~k)- I d d xJ{x)4> FP {x). (2.47) 



1 r d d k ~ 

2 J Wf a 

— IT 

The fixed point field operator is 

(b FP (x) =J2z(x~n)<P n , (2.48) 



with the coefficient function Z{x) in configuration space 



ZW = / (2.49) 



Using the fixed point field, we calculate the fixed point Laplacian in d dimensions 
for an arbitrarily placed d — 1-dimensional hyperplane as a boundary. Consider 
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zero boundary conditions <f>(x) — 0, Vx € T, where the boundary T divides the 
hyperspace into the halfspaces A inside and B outside the boundary (see Fig. || for 
d = 2). Start with the identity 




Figure 8: The line T as a boundary in d = 2. In the text we show how to calculate 
the perfect Laplacian restricted to the halfplane A. 



d d x {<j>{x) 



<f>{x)f 



(2.50) 



which is trivially fulfilled because we impose the condition <p(x) = — </>(x) on the 
field, where x denotes x mirrored at the boundary. If we plug the fixed point field 
(2.48), ( 2.49[ ) mto Ec L- ( 2.5C| ) , we get the quadratic form $ T C$ = 0, where $ is the 



lattice field vector and the elements of the symmetric matrix C are 
Cnm = I d d x [Z(x — n) + Z(x — n)] [Z(x — m) + Z(x — m)] 



(2.51) 



As it doesn't matter how we order the lattice points when arranging them in a 
vector, the lattice field vector $ can be split into the subvectors $ A that collects 
all the components Q nA living on the lattice points in halfspace A, and $ B that 
collects the components $„ B in B. The quadratic form then gets 

{<S> A ) T C A <S> A + 2(<t> B ) T C BA $ A + ($ B ) T C S $ B = 0, 



where Ca, Cba and Cb are submatrices of C 

r< _ | Ca Cab 



(2.52) 



(2.53) 



with Cba — Cjs- Keeping <I>^ fixed and minimizing Eq. (2.52) in we find the 
relation between lattice points outside and inside the boundary $ B = — Cg l CsA^ A - 
Plugging this into the perfect lattice Laplace equation , p*(n — n')$ n / = 0, the 
sum over all lattice points n' is replaced by a sum over the halfspace A 



E 



p*(n 



BAjni 



± r. 



0. 



(2.54) 
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and the perfect Laplacian p* A (nA,n' A ) restricted to the halfspace A is given by the 
term in square brackets. As the fixed point field operator is local, this exact result 
can be used to construct another parametrization of the fixed point Laplacian. 
The position and orientation of the hyperplane relative to the lattice are then the 
parameters and it only remains to calculate the couplings numerically. 



8. Conclusion 

The fixed point action has its applications not only in field theory, it can also serve as 
a powerful tool for the solution of partial differential equations, where we encounter 
non-trivial boundary conditions. We showed that the influence of boundaries on the 
fixed point action is highly local. Therefore we could provide a parametrization for 
the perfect lattice Laplace operator in d — 2 which is easy to use and gives nearly 
perfect results for any resolution. A generalization to d dimensions is trivial, as 
none of the calculations is specific to two dimensions. While the parametrization 
provided in this work can be used to approximate boundaries of any form, it is 
not the only one possible, and we proposed an alternative to find a parametrized 
fixed point Laplacian for non-trivial boundaries. The benefit from using perfect 
discretizations is enormous, as we pointed out for test problems. 
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